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ABSTRACT 


The traditional method of data presentation for heat ex- 
changer surfaces does not permit comparison of individual sur- 
face types in any simple manner. This data is most commonly 
presented in terms of heat transfer coefficients and friction 
factors referenced to the exposed area as a function of Rey- 


nolds number based on the minimum free flow area. 


Soland [5] proposed a method of surface comparison for 
plate-finned heat exchanger surfaces in which the heat trans- 
fer coefficient and friction factor is referenced to the base 
area and Reynolds number is based on the open flow area, as 


though the enhanced surfaces were not present. 


Theory and principle inherent in the Soland proposal are 
used to develop a method of surface comparison for finned tube 
surfaces as presented in (a) Kays and London [2], and (b) sev- 
eral supplemental surfaces furnished by the Trane Corporation. 
The derived comparison method is evaluated by application to 
a practical crossflow finned tubular heat exchanger design. 


Appendix III provides a method for sizing crossflow 
finned tubular heat exchangers. 


Thesis Supervisor: Professor Warren M. Rohsenow 
Title: Professor of Mechanical Engineering 
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NOMENCLATURE 


Minimum free flow dimension between tubes measured 
diagonally between adjacent rows from base to base 


as if fins were not present; see Figure 4, [FT] 


Minimum free flow dimension between tubes measured 
transversely to the direction of flow from base to 
base as if fins were not present; see Figure 4, 
Pear) 


Heat transfer area of the bare base tubular sur- 
faces; equals length times heated perimeter times 


number of tubes; [Pre] 


Exchanger minimum free flow area, defined by (la), 
[FT?] 


Constant from linearized equation for friction fac- 
tor (£) of the form: A RE ® 


Exchanger flow passage frontal area ignoring any 


enhancing surfaces; see Figure 4, [FT?] 
Ratio of fin area to total transfer area, (PT? 
Exchanger total frontal area, [Pro] 


Constant from linearized equation for Colburn mod- 
ulus (3) of the form: A RE ® 


Total heat transfer area of exchanger on one side, 
defined by 8 times total volume, [FT“] 





Jo 


Longitudinal center-to-center tube spacing, see 
Figure 4, [FT] 


Flow stream capacity rate [WC], [BTU/HR-°F] 
Constant used in the computation of friction factor 
(£) for flow over smooth tube banks; defined by 


(32) 


Constant required for the computation of radial 


fin efficiency; defined in Figure 5. 
Specific heat at constant pressure, [BTU/LBM-°F] 


Constant required for the computation of radial 


fin efficiency; defined in Figure 5. 
Hydraulic diameter, defined by (la), [FT] 
Nominal diameter, defined by (lb), [FT] 


Transverse center-to-center tube spacing; see Fig- 
ure 4, [FT] 


friction factor based on total area (Aj) i defined 
by (4a) 


friction factor based on base area [A.J defined 
by (4b) 


friction factor based on total area (An) for a 
smooth surface, defined by (31) or (32) 


32.174 [LBM/LBF] [FT/SEC] 





Mass flux based on minimum free flow area, [AJ]. 
defined by (2a), [LBM/HR-FT*] 


Mass flux based on free flow area, [AL], defined 
by (2b), [LBM/HR-FT*] 


Heat transfer coefficient based on total area, [A,. 
defined by (5a), [BTU/HR-FT*-°F] 

Heat transfer coefficient based on base area, [A]; 
defined by (5b), [BTU/HR-FT“-°F] 

Colburn j-factor based on total area, [AQ] i Gefined 
by (7a) 


Colburn j-factor based on base area, [A. li defined 
by (7b) 


Colburn j-factor based on total area, [A], for a 
smooth surface, defined by (28) 


Thermal conductivity [BTU/HR-FT-°F] 


Fin height defined by ie Pe, for rectangular fin 
sheets, see Figure 5, [FT] 


Flow length along the axis of the tubes (i.e. for 
the fluid inside the tubes); for fluid flowing 
normal to the tube banks, defined as N times b; 
[FT 


Component of fin efficiency (ney)? defined by [10] 


Number of tubes 





Nu 


Nu 


NTU 


ANcorr 


Py 


AP 


q/A 


Symbol used in reference [2] for RE 
Nusselt number; defined by (6a) 

Nusselt number; defined by (6b) 

Number of transfer units; defined by (21) 


Fin temperature effectiveness of a longitudinal 
fin of rectangular profile; defined by (9) 


Correction factor applied to Was (9) “to {perinile “the 


use Of the simplier [Tanh m2]/m2 equation for fin 
temperature effectiveness in lieu of Bessel Func- 


tion type equations. 


Total surface temperature effectiveness, defined 
by (8) 


Prandtl Number [uc /K] 

Pressure [LBF/FT?] 

Pumping Power [HP] 

Friction pressure drop [LBF/FT*] 
Heat transfer rate [BTU/HR] 

Heat Flux [BTU/HR-FT*] 


Reynolds number based on minimum free flow area, 
[AL]. defined by (3a), Labeled Np by ref [2] 
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REw Reynolds number based on free flow area, [An], de- 
fined by (3b) 

aa Radius of enhanced (finned) surface, [FT] 

ry Outside radius of base tube, [FT] 

rh Hydraulic radius, equals hydraulic diameter divided 
by four for a round tube, [FT] 

z Temperature, [°F] 

U Overall heat transfer coefficient [BTU/HR-FT*- °F} 

V Volume (FT?) 

W Mass flow rate, [LBM/HR] 

X,Y,2Z Principle dimensions of heat exchanger, [FT] 

Xp Ratio of diagonal tube spacing to tube outside 
diameter 

xX Ratio of longitudinal tube spacing to tube outside 
diameter, equals b divided by Do 

Xn Ratio of transverse tube spacing to tube outside 
diameter, equals F divided by Do 

Miscellaneous 

fo] In tube geometry, angle defined by the longitudinal 


centerline and the diagonal centerline; see Figure 





we 


7, [degrees] 


Ratio of total transfer area on one side of the 
exchanger to total volume of the exchanger, label- 
ed as a in reference [2], (pr? /Pr>} 


-0.00259, Constant used to calculate CA 


~0.00840, Constant used to calculate CA 


~0.02426, Constant used to calculate Dj 


~0.04521, Constant used to calculate D, 


Fin thickness, or average fin thickness if tapered, 
[FT] 


Exchanger effectiveness 
Viscosity [LBM/HR-FT] 


Density [LBM/FT> ] 


Subscripts 


d 


diagonal 


enhanced surface 


heat exchanger fin material 


outside 





ie? 


Subscripts (continued) 


os one side 

Ss smooth surface 

c transverse 

T Total 

2 Finned side, colder fluid flowing over finned tube 
banks 


1 Tube side; hotter fluid flowing inside smooth tubes 
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I. INTRODUCTION 


A. Purpose 


The objective of any heat exchanger comparison method is 
to enable the designer to select, from an often overwhelming 
number of enhanced surfaces, that unique surface which can be 
termed “optimum". The comparison technique should be in the 
form of a logical, easily implemented, inherently accurate 
procedure. 

It should be recognized from the onset that regardless 
of the definition of "optimum" selected, the surface labeled 
as such may not, in fact, be the best for a given practical 
Situation. Space limitations or cost restraints frequently 
restrict the designers flexibility. While a particular surf- 
ace may be "optimum" from a least volume criteria, for exam- 
ple, it may well exceed the cost limitations or require an 
unacceptably long dimension. Violation of either of these 
practical constraints could well force the selection of an 
alternate surface. 

Soland [5] proposed a method of comparison that permits 
performance comparisons of all of the Kays and London [2] 
plate finned type heat exchanger surfaces on four different 
bases? 

a. Same exchanger shape and volume 
b. Same exchanger volume and pumping power 


c. Same pumping power and NTU 





17 


d. Same volume and NTU 

Sheldon [6], (a) applied Soland's method to practical 
heat exchanger design problems; (b) evalueated the method ef- 
fectiveness; and (c) performed a comparative analysis on plate 
finned surfaces in addition to those found in Kays and London. 

The purpose of this paper is to: 

l. adapt Soland's method for application to crossflow 
finned tube type surfaces presented in (a) Kays and London; 
and, (b) several supplemental surfaces furnished by Trane 
Corporation; 

2. evaluate the method effectiveness by application to 
practical crossflow finned tubular heat exchanger design 


problems. 


B. Background 


Crossflow finned tubular heat exchangers have found app- 
lication as gas turbine plant intercoolers, aircraft engine 
and electronics coolers, air conditioning units, heat ex- 
changers, and numerous other uses. These applications most 
often involve gas-liquid service. [See Figure 1] 

The performance of a given heat exchanger in which the 
ratio of the heat transfer coefficients between the two work- 
ing fluids is greater than about three is markedly improved 
by the employemnt of enhanced surfaces. The increased sur- 


face area inherent with the employment of fins yields a net 
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improvement in heat transfer despite lower heat transfer co- 
efficients. Finned tubular heat exchangers have surface-to- 
volume ratios (8) ranging between 20-200 SO esi 

Reference [2] will be considered the primary source for 
data related to the design and testing of finned tubular sur- 
faces, while supplemental data has been provided by the Trane 
Corporation. 

Optimum design is often used synonomously with "lowest 
cost". Cost is frequently defined as (a) amortized acquisi- 
Sition cost related to the initial design and fabrication of 
the component, (b) Operating costs, defined as amortized acq- 
ulsition cost plus operating expenses and (c) cost to the 
system in terms of either weight or volume. 

Consider aircraft engine component design, for example, 
where weight and volume are valuable quantities which can be 
directly equated to dollar cost. An unnecessarily large sup- 
port component impacts either payload or endurance. 

It is certainly conceivable that a heat exchanger could 
be designed with (a) the least acquisition cost, (b) the low- 
est projected operating expenses, and yet (c) have the high- 
est weight or volume impact on the parent design and there- 
fore by unacceptable to to "impact cost". 

The estimation of dollar cost of a particular heat ex- 
changer design is difficult to estimate, primarily due to the 


proprietary nature of the required cost estimating relations. 
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Companies that design and manufacture heat exchangers are of- 
ten quite reluctant to release this data. 

A method of dollar cost comparison will be proposed, but 
due to lack of data, will not be fully investigated. 

The term "cost", therefore, as used in this paper, will 
most generally refer to "impact cost" and will serve as a 
means of relative ranking of proposed designs specifically by 


required volume, pumping power and NTU comparison. 


C. Surface Testing and Data Presentation 


In order to provide the required friction and heat trans- 
fer data for a specified finned tubular surface, experimental 
testing is required. 

The experimental procedure employed by Kays and London 
is clearly outlined in reference [2]. 

Regardless of the technique, the results include fanning 
friction factor f, Prandtl number, Reynolds number, and Col- 
burn j-factor. Figure 2 is an example showing the tradition- 


al method of data presentation. 


D. The Task of the Designer 


Assuming that the designer of a heat exchanger has de- 
fined such constraints as: 

(a) allowable pressure loss, 

(b) temperature change, 


(c) amount of heat to be transferred, 





eu 


FIGURE 2: Sample data plot showing traditional method of data 
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(d) volume and weight restrictions, 

(e) fouling and corrosion considerations 

(£) materials 
he is now ready to select the heat transfer surface(s) to be 
employed. 

Having compiled a great number of candidate surfaces 
from (a) existing data, or (b) having designed and tested his 
own surface, he is now faced with the task of comparing these 
surfaces to determine which should be selected for the specif- 
ic application under consideration. 

Figure 3 presents two examples from Ref [2] showing the 
traditional method of presenting f and j verses RE data . Note 
that in each case, Surface B has a superior heat transfer 
characteristic at a given Reynolds number; however, Surface A 
has a lower (and therefore, superior) friction factor at that 
same Reynolds number. Which surface is "optimum"? 

There appears to be no obvious answer to this question 
and this fact has prompted several methods of surface compar- 
ison to be developed. 

It is intended that the method presented herein will pro- 


vide a solution to this dilemma. 





0645" 
i!49' a fe es: 0.0!0" 
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Surfaces of Ref [2] Figure 97 
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FIGURE 3: Two examples from REF [2] showing the traditional 
method of presenting f and j verses RE surface data, 
and showing the difficulty in comparing surfaces A and 


B in each case. 
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II. PROPOSED COMPARISON TECHNIQUE 


A. Derivation of Basis of Comparison 


Soland, in reference [5], provides a detailed derivation 
of the proposed surface comparison technique used to analyze 
the performance of crossflow plate-finned heat exchangers. 

Theory and principle inherent with the Soland method are 
used to derive a comparison technique to analyze the perfor- 
mance of crossflow finned tubular heat exchangers. 

Comparison of the performance of various finned tubular 
surfaces assumes the following quantities are held constant: 

A) W, flow rate 


B) T hot fluid inlet temperature 


h,in’ 


C) T_.., cold gas inlet temperature 
c,in 


Other important assumptions are: 

E) The heat transfer resistance of the tube 
walls is negligible. 

F) The controlling heat transfer resistance 
is assumed to be on the finned side of the 
heat exchanger. 

Traditionally, data is presented in terms of heat trans- 
fer Colburn modulus, j, and friction factor, f, based on the 


total exposed area, A as a function of Reynolds number, RE, 


[ue 
based on minimum free flow area, Aas and a hydraullic diameter, 
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Dis of the flow passage. 
The proposed comparison method converts these j and f mag- 


nitudes to new quantities In and f,, based on the heat transfer 


N 
area of the bare base tube surface, Aye and a Reynolds number, 


RE which is referenced to the open flow pasSage area, A as 


N’ FF! 
though the fins were not present. The effect of the fins is 
accounted for as an increased heat flux and hence larger h for 
for the bare base tube surface area. In order to incorporate 
the effect of the fins into Jy the metal conductivity must be 
specified. 

Table I shows the proposed new definitions and compares 
them with the analogous definition used by Kays and London [2]. 
Two cases must be considered as illustrated in Figure 4. 

In order to convert data presented in the format of ref- 
erence [2] to the new basis, various ratios are derived from 
basic definitions, equations (1) through (10), and Figure 4. 
These ratios are presented in Table II. 

Two assumptions and an explanation of the calculation of 
radial fin efficiency are required in order to solve for the 


proper fin efficiency, Nye needed in the conversion relations: 


~ > fj 


1. The results presented herein will assume the gas 
flowing over the finned tube banks is air whose 
thermodynamic properties will be evaluated at 90°F. 

2. Some of the results presented herein specify copper 
fins [K=222 BTU/HR-FT-°F] and others specify Alumin- 
um fins [K=117.6 BTU/HR-FT-°F]. 
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CASE I: Minimum Free Flow Area Occurs 
In Spaces Transverse To Flow 
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CASE II: Minimum Free Flow Area Occurs 
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FIGURE 4: The Two Cases of Minimum Free Flow Area Occurance 


and the Associated Geometric Relations of Each CaSe. 
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3. Figure 5 explains a proposed generalized procedure for 


estimating 
provides a 
in lieu of 
tion type 


Figure 6 shows 


radial fin efficiency. This procedure 
Simplifed method of efficiency calculation 
uSing more technically correct Bessel func- 


relations. 


an example of data presented on both basis. 


It should be noted that additional curves of Jn verses REY 
would result for other magnitudes of fin thermal conductivity. 
For any heat exchanger, the fluid pumping power per unit 


| volume on one side 18S given by: 


= S (17) 





From equation 4b: 


2 
AP, _ tN 4 N b Gy 
Diy 20 Do 
From equation 2b: 
W = Gy Ap 
3 3 
2 i SE FRE b 
= lala: rao" Ma (18) 
os 





Radial Fin of Rectangular 
Profile 













8Neorrection” 


where: 


c,° ~.00259 - .0084 in 


elo" elo 


o,° ~.02426 = .4521 In 


{Ref 3 page 264] 








Compute fin efficiency assuming a 
Longitudinal fin of Rectangular 
Profile: 


. tanh me 
ay m 


What geometric shape describes the 
fin being employed on the tube 
sutface? 


Tapered (Hyperbolic) Rad- Rectangular Profile fin 
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[see below for £ defn) 
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FIGURE 5: Computing fin efficiency of a Radial fin 
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FIGURE 6: Comparison of Colburn j factors (j and Jig) and com- 
parison of friction factors (f and £.)) for the Sur- 


face of Figure 94 [Ref 2]. 
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For the same fluid at the same temperature level, tu and 0 


are constant: 


5 
p £, RE, A, b 
a cl) 
Vos Be Av 
N os 


Figure 7 shows the geometry involved in the calculation 


of [A,, bI/V.. for both case I and case II. 


The heat transfer for any heat exchanger is given by: 


q=e(T i. - To inlW c, (20) 


For any given flow arrangement, ¢€ can be related to NTU 


by either a closed form equation or graphical curve similar to 


that shown in figure 8, where 


Ah 

NTU = ——— (21) 
Wc 
Pp 


ra 


The relationship between e« and NTU is always monotonical- 
ly increasing. In other words, an increase in A hy causes an 
increase in NTU which means e€ and hence q are greater, given 
that the fluid properties and flow rate are held constant. 


From equations (21) and (7b): 
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Case I: Minimum Free flow area Transverse to Flow Direction 


Volume on one Side: 


T p* 
40 fe) = 


L dimension into paper along 
flow length of tubes. 


[X,D,} [X,D JL = 


IN b" Volume: V/LL/A 


Ap b= [XpD, - Dit xy De 
Ratio: 
Ap b : [X,, ~ 1)x, 
>IT 
V XpX 4 


OS 





Case II: Minimum Free Flow area in the Diagonals 





Total Volume on one Side: 


. 2 
V = [D2 X 


: TE 
os ~~ Zo! 


2 


A, b" Volume: * 


A, b 2X, Dy [X, eel). leg 
Ratio: 


Ap b 41X, - 1]X, @ =D 
Vos _— - 7 L/ 


FIGURE 7: Calculation of (A, bI/Vo5 for both case I and II. 
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A, h eer Ge, C 
ntu _ “pb "Nn _ *b In Sn Sp (22) 


V V wc 23 
os os p ae Ve W CS 


From equations (3b) and (lb): 








wry _ 4 Ap B Jy REY U 
os pD. pet wv 
Nery os 
Op a 
mn SG 4 RES Aa b 
23 
Os Be Dy os 


For the same fluid at the same temperature level, Cy Uy, 


and P are constant: 








RE A, b 
5 25) 
N 


And for the flow rate, W, also constant: 


j 
ee N "NF (26) 
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NTU 


FIGURE 8: Typical Graphical Relationship between heat ex- 
changer effectiveness and NTU. In general it is 
possible to express: 


/C 


e = O(NTU, C flow arrangement] 


min’ ~max’ 














Si 


Recognizing that the f£ and j verses RE plots for flow 
over finned tube banks from reference [2] can be closely ap- 
proximated by a straight line on a log-log plot, implies that 


f and j can be represented by an equation of the form: 


fs 


. 
z 


(e) °f 


J(e) is bal 

Table I of Appendix III lists the linearized equations of 
f and j for the surfaces represented by figures 92 -101 of 
reference [2]. This simplification allows for ease in pro- 
gramming the derived surface comparison equations as demon- 
strated in the procedure outlined in Appendix I. 

In any event, the performance parameters may be calculated 
from enhanced surface data once the data is presented in the 
This presentation may 


form fA verses RE,, and Jn verses RE 


N N° 
be (a) graphical due to the complex relations presented by f 

and j} verses RE, as suggested by Soland or (b) analytical, as 
outlined in Appendix I, due in large part to the linear nature 


of f and j verses RE on a log-log plot. 


B. Method of Surface Comparison 


Having expressed the enhanced surface data in the form fy 


verses RE, and Jn verses RE it is a relatively simple matter 


N N’ 
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to calculate the performance parameters of equations (19) and 
(26) and plot them. Figure 9 1s an example of such a plot. 

The performance parameters of two hypothetical surfaces, label- 
ed 1 and 2, have been plotted in order to demonstrate how a 
determination of heat exchanger relative performance may be 
made. 

Four different comparisons are immediately available from 
figure 9 and are indicated by points a, b, c, and don surface 
2. Point o on surface 1 represents the reference heat exchanger 

design to which each of the four points on surface 2 will be 


compared. 


Case A: Same heat exchanger shape and volume (Li=Ly: Viger 
A, =A, ). Because W and A, are fixed: 
a re) 


{D,, ] 
RE. = RE es (27) 


a 
No Ne [D,, 

O 
The results of this comparison are easily obtained as 
the ratios of ordinate values [NTU/V. 2] and abscissa 


values R/V ae of the two surfaces. 


Case B: Same heat exchanger volume and pumping power (Vi=Vo. 


PH Po) > Point b is located on a vertical line through 


point o since the pumping power per unit volume is 
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FIGURE 9: 


3 
fy RE AL b 


os Diy a7 


belly 
R 


Performance Parameter Curves For Two 
Hypothetical Surfaces Showing Points 


Used in Sample Comparisons. 
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Case B: (continued) 


identical in both exchangers. The NTU ratio of the 
two heat exchangers is obtained as the ratio of the 


Ordinate values from Figure 9. 


Case C: Same pumping power and number of transfer units. 
(P =P oe NTU =NTU.) Point c is located on a line hav- 
ing slope equal to unity and through point o since 
both NTU and P are constant and each axis is inverse- 
ly proportional to volume. The ratio of the volume 
required using surface 2 to the volume required using 
surface 1 is the ratio of either ordinates or abscis- 
sas at points o and c, respectively. 
Since material costs are usually a function of volume 
Or material weight, and thus heat transfer surface 
area, it is through Case C that a relative cost com- 


parison may be attempted: 


[Total Cost], [$/FT“], [8], [Vy 


aay 4 
DI La? uc... Gao 
{Total Cost], [S/FT FT (8), [Ve], 


Note that as long as surface 2 lies above surface l, 
the result will be a smaller volume required to do 


the same "job". {see Figure 9] 


Case D: Same volume and number of transfer units. (Vi=Voe 


NTU J=NTU or qg=4,) > This case compares the pumping 
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Case D: (continued) 
power required by surface 2 to that required by sur- 
face 1 for the case when both heat exchangers yield 
the same overall heat transfer performance. The com- 
parison is made by taking the ratio of the abscissas 
(points o and d@) in Figure 9, since a horizontal line 
on the performance plot has a constant value of 
NTU/V.,° The surface 2 Reynolds number for this case 
is the smallest of all four cases and consequently, 
the flow area is the largest. Since the volume is the 


same as the surface 1 exchanger, the length must de- 


crease. 


C. Summary 


Using the equations derived herein it is possible to con- 
struct a performance parameter plot of NTU/V.. verses B/ Vs 
(see Figure 9), from which it is possible to obtain useful 
performance comparisons between two heat exchangers for four 


different criteria. 
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III. COMPARISON OF FINNED TUBE SURFACES 


A. General 
Using the procedure derived in Chapter II and outlined in 
Appendix I, all the finned tube surfaces were plotted as per- 


formance curves of: 


: 3 

3,, RE b f.. RE. A, b 

» N “P verses A ee. 
Dy Oe Dy Vos 


These performance curves assume constant W, T. Bane and 
¢ 


T » and are equivalent to NTU/V. verses P/V. 


h,in Ss. 

Figure 10 shows the resulting performance parameter plot 
presenting only the highest performance curve at each nominal 
diameter [Dy1- This study considered sixteen different sur- 
faces having eleven different nominal diameters. Calculated 
performance results are presented for each surface considered 


in Appendix II. 


B. Comparison to Smooth Surface 


All of the sixteen finned tubular surfaces were then com- 
pared to a smooth surface having no enhancements (fins) and a 


nominal diameter identical to that of the enhanced surface. 
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PERFORMANCE PARAMETER CURVE FOR FINNED TUBULAR SURFACES 


FIGURE 10 





10° 10°? ro24 1022 1023 10 


3 
Pp £..RE Ap b 
os pf V 
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Note that subscript e will be used to denote the enhanced sur- 


face and the subscript s will be used to denote the smooth 





Ss es ee ee | 
° ’ 
SHO h ce Js Sy 25 e calculations were completed 
° p2 v | 
: N “os 
b. 
using eh skkas [7] correlati for forced-convection, 
® 
turbuldh f£ rer staggered tube*fbanks with the appropriate 
” fue 
oe ameter [DJ - 
 ) 
: 
IQ justify the use oad the, above relation to 
»x 





predict ShoRth™ t be bank Colburn Mddulus [j_], equation (28) 
@ e 

was usge toto ct the actual Reat ¢ransfer data for flow 

over barg ited nks as presented Hf, reference [2]. The re- 


sults re! shown n figure ll. Makan mum deviation between the 
v' e oe 


exper imentai val e of a (from reference 2) and the calculated 
A 


value of +4 tfror equation (28) is 47.58. , 
" \ ¢ 


Note: that; jL may be presented] for each individual enhan- 
r .4 a 


Poe =f 
(% co” . 
ced surface, in’ the form: | 
: an a 


n 3 ; 
. eh By (3) * ' 
Ig * A5(s) REy ae) 


For a smooth surface h=hy- G=G 


a 
C rn 


we RESRES 
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Equation (26) may then be written for the enhanced sur- 


face as "if smooth”: 


NTU. 2 Js REw Ay b 
—o oa 
os D V 
N os 
f_ RED A, b | 

Smooth surface ‘s N F were completed using a corre- 

4 

Diy Vos 


lation for friction factor (f£.) derived from first princivles 
using the Kays and London [2] data for turbulent flow over 


staggered bare tube banks as the data base: 


Case I: 
— C 
£ = 0.7184 REDO'7*49! py ~339° 793 (x 3 7 (31) 
5 * L 
where: C,, = 0.828 - 16.596 re 04277 
Case II: 
— C 
£ = 1.0000 RE~?° 1892 fy ~94 9-793 fx 3} P (32) 
5 D L 
~0.7284 


where: Cy = 130.10 RE 


Reference to Case I and Case II differentiates between 


minimum free flow area occuring transversely (I) or diagonally 
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(II) to the flow direction. (see Figure 4). 

In order to justify the use of equations (31) and (32) to 
predict smooth tube bank friction factor (f.), these relations 
were used to predict the actual experimental friction factor 
data for flow over bare tube banks as presented in reference 
[2]. The results are shown in Figure 12. Maximum deviation 
between the experimental value of f. (from reference 2) and 
the calculated value from either (31) or (32) is =3.03%. 

The friction factor, for each individual enhanced surface 


may be presented in the form: 


B 
RE f(s) 


fs = Ae(s) REy (33) 


Equation (19) may then be written for the enhanced sur- 
face as “if smooth": 


3 
£ RE. A, b 
aS N F (34) 


Of the sixteen surfaces investigated, five have a nominal 
diameter [D,]=1.60 inches. Since this offers an excellent 
means of comparison, cases A,B,C, and D where calculated for 


these five surfaces: 





Ge Gaus e 
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‘ o 9 g 2 Bo, °° 2 
SS Oe ae or ae A ce cn a a 7 00°0 
co is] a9 

S0°0+ 
OT°O+ 


*{Z) Jeu uT paejusssad se (>3) AOQIeJ uoTtWOTAZ jo 
OnTeA peuTtejzqo ATTejueuTzedxse butjZoOtpead ut syueq 
eqn3 perebbeys zdsao MOT AUeTNqQainy *uOT39D90AU09 
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Figure 13: 


Figure 14: 


Figure cS: 


Figure 16: 
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Shows the ratio of ordinates [NTU /NTU _] 
verses RE, for Case A. [Same heat exchanger 
shape and volume] Since both the smooth 
and enhanced surfaces have identical nomin- 
al diameters, the nominal Reynolds number 


is also identical for each surface. [see 


equation 27; results read from Appendix ITI] 


Shows the ratio of abscissas Craze) verses 
REw for Case A. Since both the smooth and 
enhanced surfaces have identical nominal 
diameters, the nominal Reynolds number is 


also identical for each surface. [see 


equation 27;results read from Appendix II] 


Shows the ratio of ordinates [NTU /NTU _] 


verses REw for Case B. [{P=constant and 


=constant ] 


Shows the ratio of abscissas or ordinates 


to produce Wie Vel verses RE Recall that 


N° 
each axis of Figure 10 is inversely pro- 


portional to volume. This is Case C. 
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Figure 17: Shows the ratio of abscissas [P_/P 1 
verses REw for Case D. [NTU=constant, and 
V=constant |] 


Higher ratios are preferred for Figures 13,15,16, and 17 
while a lower ratio is most desirable in Figure 14. 
As a matter of secondary interest, of the five surfaces 


with identical nominal diameters, there are three discrete val- 


ues of the ratio [A,, bI/Vo5 


Surface [AE bI/Vo5 

32.66; 93 1.4282 

96 & 97TA ip Bo9 
100 1.4014 


Comparison of many surfaces having identical nominal dia- 
meters to a single common smooth surface having the same nom- 
inal diameter permits the best relative comparison between en- 


hanced surfaces. 








55 


15,000 100 


10,000 
9,000 
8,000 


7,000 
6,000 


3,000 


4,000 


lm 


3,000 


rg 


2,000 


1,500 


1,000 
900 
800} 


3 2 3 ey 5 Gane? 662 


10 10. 


REw (smooth surface) 


FIGURE 17: Performance comparison results for CASE D [i.e. NTU 
and V are constant] Di =~ 1.60 inches for both the 
enhanced and plain surfaces. 
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IV. HEAT EXCHANGER DESIGN PROBLEM 


A. Description 


Appendix III is a procedure which can be used for sizing 


crossflow finned tubular heat exchangers. 
taken from reference [2], 


heat exchanger size for the given conditions. 


The following data, 


is used to determine the required 


The design is that of an intercooler for a 5000 SHP gas 


turbine plant. 


DESIGN CONSTRAINT 


Surface 


an 


6 


SIDE I 


Smooth Tube 
0.0625 FT 
N/A 
5.9291 Fr¢/rr3 
N/A 
N/A 
400,000 LB/HR 
60° F 
82° F 
1.554 PSI 
2.360 LB/HR-FT 
1.000 BTU/LB-°F 


62.38 LB/FT? 


SIDE ITI 


Finned Tube [98C] 
0.0445 FT 
0.0010 FT 
61.9 FT°/FT° 
0.8350 
0.3445 IN 
200,667 LB/HR 
260° F 
80° F 
0.3060 PSI 
0.0482 LB/HR-FT 


0.2436 BTU/LB-°F 


0.1755 LB/FT> 
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DESIGN CONSTRAINT SIDE I SIDE II 
Ey 6.8000 0.7000 
Ka N/A 222 BTU/HR-FT-¥F 
F 2e725- IN N/A 
b 1.750 IN N/A 
A, 0.4418 IN? Calculate 


Using the design procedure outlined in Appendix III, the 


principle dimensions are: 


X 6.43 FT 
Y°=e26.11 FT 


Z 


S63 21 


Total Volume = 105.46 FT” 


This is a possible heat exchanger design that satifies 


the given constraints. 


B. Attempt to Minimize Total Volume Be 


If minimization of total volume is a required design con- 
straint, the engineer has several options the limits of which 
might be bounded by: 


(A) Design a heat exchanger to the given specifications 
using each known surface for which f and j data are 
available. In this study sixteen surfaces are con- 


Sidered. 


(B) Noting the marked superiority of reference [2] figure 
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(B) continued 


94, apply the method developed in Chapter II to es- 
timate total volume savings relative to the base or 
reference heat exchanger designed in part I. 
Option (A) will be eliminated as being impractical due to 
the required effort. 
Option (B) will be shown to be the most direct and highly 
accurate approach. | 
0 
The f and j verses RE data for reference [2] figure 98C 
and 94 are presented in Figure 18. The problem presented cor- 
responds to Case C described in Figure 16. 
Using the heat exchanger design procedure outlined in 
Appendix III and replacing the finned tube side with the sur- 


face of figure 94 [Ref 2]: 


X = 1.50 FT 


ef LO. 89 ET 


1... ORs 


Total Volume = 28.68 rr? 


Z 


A 73% reduction in volume would result by employing the 


surface of figure 94 [Ref 2}. 


In order to use the comparison method developed in Chap- 
ter II to quantitatively predict the volume savings, a perfor 
mance parameter plot is constructed for each surface consider 
ed. Figure 19 presents the only two surfaces being considerec 


in this example. 
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FIGURE 18: £ and j verses RE data for the two surfaces con- 


sidered in the design exampie. 
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As outlined in Chapter II for Case C: 


(1) From the design point of the reference heat 
exchanger (Figure 98C) graphically or analyt- 
ically construct a line of slope l. 


(2) The intersection of this line with the new 
surface (Figure 94) will yield the performance 
parameters at the design point of the new 


heat exchanger. 


(3) The predicted volume reduction is found from 
the ratio of either the ordinates or abscissas 


of the two surfaces. 


The performance parameters of the reference heat exchang- 


er design are: [Reference 2 figure 98C] 


= = 401.20 [1/IN7] 

Os 
—_ = 2.93 x 10°* [1/1N*}] 
OS 


The performance parameters of the heat exchanger that 


would result from using the surface described by reference [2] 


figure 94: 
are = 1634.68 [1/IN*] 
OS 
= 1.19 X 107? [1/1N‘] 


OS 
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In other words, a predicted volume reduction of 75% is 


anticipated using the comparison method developed herein. 


C. summary 


Two surfaces have been compared from a minimization of 
total volume point-of-view. It has been shown that the per- 
formance parameters developed in Chapter II allow us to com- 
pare any number of surfaces without going through a complete 
heat exchanger design calculation. Where there may well be 
hundreds of different surfaces to be considered, the value of 


such a comparison technique cannot be overestimated. 
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V. CONCLUSIONS 


1. The comparison method of Soland [5] was used as a basis 
in developing a comparison method for finned tube heat ex- 
changer surfaces. 
2. Heat Exchanger Performance can be compared on four 
different bases: 

a. Same shape and volume of heat exchanger 

b. Same exchanger volume and pumping power 

c. Same pumping power and NTU 

ad. Same volume and NTU 
3. All of Kays and London [2] finned tube surfaces and 
several surfaces supplied by the Trane Corporation were com- 
pared on the above basis. 
4. The method developed herein was applied to a practical 
heat exchanger design problem and found to accurately predict 


relative volume savings between different surfaces. 
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APPENDIX I 
A. Purpose 

A frequent criticism of academic solutions to 
"optimization" of design is: "So what, .. . how does this 
help the engineer in industry?" [8] 

It is hoped that this section will give the "engineer 
in industry" the needed background and technique to effective- 
ly compare finned tubular heat exchanger surfaces. 

A step-by~step flow diagram using the performance 
comparison equations derived in the text of the thesis is 
presented in Figure A-l. This procedure is easily programmed 
for use on either a (1) standard computer or (2) programmable 
calculator. 

In keeping with the spirit of the criticism of 
Reference [8], a program for use on a personal programmable 
calculator will be presented herein. It is intended that 
this simplified, easy-to-use approach will demonstrate that 
the “engineer in industry" needn't resort to often costly 
computer time/programs in order to effectively "optimize" 


a design. 


B. Proposed Program(s) for Calculating Surface Performance 
Parameters on a Programmable* Calculator 


The following programs follow the general outline 
Of Figure A-l: 
Figure A-II: For calculating Surface Performance 


Parameters for the case in which the 
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CALCULATION OF FINNED TUBULAR JIEAT EXCHANGER SURFACE PERFORMANCE PARAMLTERS 






INPUT 


Linearized equations of the 
form A[RE]} for £ and } for 
both the enhanced and “if 
smooth” surface: 







A 





ste)! Bs (e) 


Recs); Pe(sy? 4¥(s)* 35 (s) 


Geometrical Factors: 
Xi Xn} Dot Dopi Sepa ce: 


8 ; [A Lame) i 


Other: Air properties, fin 


Acca)? Bee)? 







fin 





material and shape, 
Revnolds Wumber range 







Was 







DOES THE MINIMUM FREE FLOW 
AREA OCCUR IN A DIAGONAL? 










Calculate the Calculate the 
Nominal Diameter: 


4D_X, [X,- 2] 








Nominal Diameter: 
_ 8D LX, (X_-1) 


T 














Compute the ratio of friction 







(ft) 161x,-117D,, factor based on base Area {A,) 
f - ee to the friction factor based 


on total area [A,) 











Compute the ratio of Flow 
[Ap } tims equivalent flow length 
{b] to the Volume on one side 


















A,,b 4{X)-1] X 


Apb 2(X.-1) Xy i 


va.» X 7X70. 25m 


os T 










os X,*sina-0.251 
[Vog! of a heat exchanger 






Compute ratio of Reynolds Num- 
ber based on minimum free flow 


are (AA) to Reynolds Number 





based on free flow area {A,] 


Compute Power per unit Volume 
on one side of a heat exchan-~ 
ger using the specified ENHANCED 


surface. B 
f (e) 
. of = RE 
Note (e) Ae (e) 





Compute Power per unit Volume 
on one side of a heat exchan- 
ger using the specified surface 
as “if SMOOTH" 





B 
Ten, x £(s) 
Note: f(.)= Aggy RE 





RE= RE,, for smooth 
surface 


FIGURE A-l: Performance Parameters 





Radial Fin of Rec- 
tangular Profile 


*"corcection” 


Cc, = % delay _)) 
atherer 


¢,° =.00259 = .6004 in 


oo.” =,02426 = .4521 In 


(nef 3 page 164) 
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ompute fin effic- 









iency assuming a 
Longitudinal fin of 








What geometric shap 





describes the fin 
being employed on 





the tube surface? 










Tapered (hyperbolic 
Radial fin 










a = 
‘neorrection 


0 


(Pat @ page 3-116} 






s-corrected_ 


correction 













Compute Surface n: 


A>: 
4in 
n el- eet gees! 














Add selected incre- 
ment to RE. Test to 


sec if range of 
Gata is exceeded? 






The fin material for a giv- 
en surface is known. Assume 
a gas and a operating temper- 
ature in order to calculate 
fin efficiency for comparison 


purposes. 





Rectangular Profile 





fin sheet 





es 
8% correction 







supe 28 Re siel Tine 
Pectanguler Profile 
using fr istaggered bar.&) 


-$6 
ete 24) Ka ue 
t. e 2 






La r -fr 
(Pef 4 page 3-116) 






Compute number of NTU per unit 
Volume on one side of a heat 
exchanger using the specified 
ENHANCED surfacc. 
Note: Jy7 65n Dy) /P oy) and 

: B. 

= A. RE Jfe 

J(c) Aj (e) 
Compute number of NTU per unit 
Volume on one side of a heat 
exchanger using the specified 
surface as “if SMOOTH”. 


B 
ee a j(s) 
Note: j/5)* AS (5) RE 


YES END 





Figure A-II-A: 


Figure A-II-B: 


Bogune A-TII: 


Figure A-III-A: 


Eagure A= Ero: 


C. Operation 


step 
i. 


Ze 
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minimum free flow area occurs transver- 
sely to the flow direction. 

Definition of geometric constraints/ 
constants to be stored in memory. 


Sample output for transverse case. 


For calculating Surface Performance 
Parameters for the case in which the 
minimum free flow area occurs in the 
diagonals. [see Figure 4] 
Definition of geometric constraints/ 
constants to be stored in memory. 


Sample output for diagonal case. 


Description 


Select and Load Program 


Store Constants as outlined in either 
Figure A-II-A or A-III-A. 


Press Key A 
Compile Output 


Stop Program when range of RE data is 


exceeded. 
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PROGRAM* FOR CALCULATING SURFACE PERFORMANCE PARAMETERS 

FOR A CROSSFLOW FINNED TUBULAR HEAT EXCHANGER IN WHICH 

THE MINIMUM FREE FLOW AREA OCCURS TRANSVERSLY TO THE FLOW 
DIRECTION 


FIGURE A-T1 
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wie Ree Sele 678 LN ae 134 1 01 198 : 24 
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FIGURE A-II-A 


GEOMETRIC CONSTRAINTS/CONSTANTS TO BE STORED IN MEMORY 
—  — 


[Transverse Case] 


neo ory DESCRIPTION ara gry DESCRIPTION 
leeve Blenk. Thie regieter ie 

e an weed in the program to etore the || 16 Ri/M Ratio of fin erea to total 
welve of RE being used in eech eree. 
iteretion. 

Veom linesrized equetion for 
2 Ber : fsiotion factor (f) of the en- 11 L Pin length [inches] 
~ fenced eurfece of the form: 
ange ® 
Froe linearized equation for Property constent defined by: 
fciction factor (f) of the en- q a/) 
-! Pete) = ‘ [C vl/ tk. Pe) CxO) 
panced eurface of the fora: P a Ff 
A RE 
From linearized equation for From lineerized equation for 
3 Mey ; friction fector (f) of the if 13 Ay ie) Colburn modulus (4) of the en- 
° gmooth eurfece of the form: hanced surface of the formar 
E ; an ® ae ene * 
Prom linesrized equation for From linearized equation for 
4 %, feiction factor (f) of the ig 14 By ce) Colburn modulus (j) of the en- 
- &(e) gmooth eurfece of the form: hanced surface of the form: 
ans ® are © 
Specified hydreullic diameter Bpecified fin thickness or 
: 3 By, of the flow pecsage [inchee] 15 é ‘everege thickness, if tapered. 
(inches) 
from linearized equation for 7 
Colburn modulus (4) of the "if 
bd 4, a ste) smooth® surface of the form: 
Coastents derived from graph on De ee: x he ® 
pege 246 of en . They . 
4 meteugse, () From linearized equacion for 
are used in the computetion of e 
Colburn Modulus (j) of the °if 
? 6; eorrection factor applied to 17 PY 
: $(s) smooth” gurfece of the form: 
. a im erder to permit the use a B 
(ob A RE 
g of the elaglier Tenh m2 equation 
C . Input one: 

6 é, fox fia teaperature ef fective- 18 Pin ID 0.000 if fin te tapered 
meee in lieu of Beesel function 71.000 if fin is @ radial fin 
type equations. of rectangular profile 

; Ratio of outside tube redius 

9 a, 19 B to outside fin radius. 


Any Cesired number to permit the 
a feta 2D user to identify the eurfece be- 
ing processed. 


Longitudinal tube pitch. Ratio 
o x, of longitudinal tube spacing to 
tube outside diameter. 


Opecified heet trenefer area per 


8 8 2 3 g of transverse tube spacing to 
totel vol wetio. fin’ /in 
: — - ( / ; } *% tube outside diameter. 
Desired increment of REL: A 
ified ¢ outside Ciame . 
Cc o Spec ° ube . ¢ -— z Increment value of 10,000 is recommnedcd. 
finches) 


: Traneverse tube pitch. Ratio 








fe 


4 i 
, 
| - 


leit 





4 


et) 
Til 








EEGURESA=1 i=B 
SAMPLE OUTPUT 


{Transverse Case] 


FIGURE ID jeep 100.0000 


0.0000 

OF Z23 

=J2073 

0.3514 

~0.0896 

0.1430 

0.0054 

-0.0084 

0.4098 

ac spat 

1.6402 

£ y/Ft ee 14. 4115 
10.2285 

2.1542 

(A, by/V, 2 1.4014 
977.6564 


Print out of Regis-~- 


ters 0-9 and A-I 


‘= 


ye) 


BE. / RE 


HmMmHoOWPr WHOA UW BWNEFH OC 


0.8390 
0.3240 
0.0001 
0.1650 
~0.3986 
0.0130 
0.5830 
-0.3960 
-1.0000 
0.3829 


Print out of Regis- 


ters 10-19 and A-I 


wo MA IHAMNPWNF OS 


Value of RE), being, 
processed 
977.6564 Value of RE being 
10000.0000 processed 
8.186009734+10 <=————- _—&P Be 
2 981639001410 
0.9328 &————-— a 
NTU./Vos — nn 598.0269 
79.1432 —————_ NTU Pig 


OS 


OSs 
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FIGURE A-III 

* 
PROGRAM* FOR CALCULATING SURFACE PERFORMANCE PARAMETERS FOR A 
CROSSFLOW FINNED TUBULAR HEAT EXCHANGER IN WHICH THE MINIMUM 


FREE FLOW AREA OCCURS IN THE DIAGONALS 





@@1 sL6Li per 3 ci 
-- ae ey 113 Z 03 165 ENT? a 
, a 114 vy 3} 170 Li! 22 
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an : — 119 RCL 76 oe | 175 A -55 
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a: ee 126 yr 31 164 ! 81 
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819 =e 138 RCL? € 3 1e€ FCLD 9-2 63 
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we: : s 158 : -24 206 = RCLC 3€ ey 
+ ~~ 151 2 1€-5} Por FCG 3€ 12 
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ae ee 153 <6 205 22 1€-51 
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ae ee ee 15d -RCLO 36 14 214 RCLE 4 Ce 
os a 155 CHE 22 215 ¥ -25 
earn: ei 168 ev 73 Big © iaebe clea 
aa : e 11 S 4 oie 2 €3 
re af je LST ie: Bye : =24 
ee pre ae 163 CLE 3€ 14 219 RCLE 7 15 
a53RCLE OE 0S io ae ears re 
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*For Hewlett-Packard 67/97 Programmable Calculator 
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FIGURE A-III-A 


GEOMETRIC CONSTRAINTS/CONSTANTS TO BE STORED IN MEMORY 


BEG arty 
6 % 
i Aete) 
~ Pete) 
os Asie) 
: ‘ Pete) 
: Pan 


« é, 
? ay 
| | ay 
9 4, 
A Date ID 
5 . 
¢ d, 


[Diagonad Case] 


DESCRIPTION 


Specified diagonal tube pitch. 
Ratio of diagonel tube epacing 
to the tube outeide diameter. 


Prom linearized equetion for 

friction factor (f) of the en- 

hanced eurtece of the form: 
arn 5 

From lineerized equetion for 

friction fector (f) of the en- 

henced surface of the form: 
are ® 

Prom lineerized equation for 

friction factor (f£) of the it 

emocth surface of the form: 
an 


Prom lineerized equation for 

friction fector (f) of the if 

gmooth eurfece of the form: 
ane? 

Specified hydraullic diemeter 

of the flow paseage (inches) 


Conetante derived from graph on 
pege 264 of Reference (}}. They 
are used in the computetion of a 
oorrection factor applied to 
Moon in order to permit the use 


of the simplier Tenh nl equetion 
ery at 


for fin temperature effeotive- 
pees in lieu of Beesel function 

° ¢t 
type equetions. 


Any desired number to permit the 
veer to identify the seurfece be- 
ing processed. 


Specifiod heat trenefer area per 
total volume ratio. (in? sin?) 


Gpecified tube outeide diameter. 
(inchee) 


REG 
10 
11 


12 


13 


14 
15 
1é 
17 
18 


19 


A/Ry 


est 


Ay (e) 


By ie) 


Aye) 


85 (6) 


Pin Id 


* Increment 


DESCRIPTION 


Ratio of fin aree to totel 
aree. 


Ratio of fin length to fin 
thickneaa. (1//inch) 


Property constant defined by: 


a/t 
(Co wls (RL Pe tw) 


Prom linoarized equation for 

Colburn modulus (j) of the en- 

hanced eurfece of the form: 
are ® 

Prom lineerized equetion for 

Colburn modulus (j} of the en- 

surface of the form: 

Spee 
Traneverse tube pitch. Ratio 
of transverse tube specing to 


tube outside diameter. 


hanced 


From linearized equation for 
Colburn wpodulus (j) of the *°if 
smooth® eurfece of the form: 
A RE > 
Prom lineerized equetion for 
Colburn Modulus (4) of the °if 
emooth® surface of the form: 
are ® 
Input one: 
0.000 if 


-1.000 if 
e of 


fin ie tapered 
fin ie a redial fin 


rectengular profile 


Ratio of outside tube radius 
to outside fin radiuve. 


Longitudinal tube pitch. Ratio 
of longitudinal tube spacing to 
tube outaide diameter. 


Angle defined by iongitudinel 
centerline end diagonei cen- 
terline: (degreos] 

Oesired increment of RE, A 
velue of 10,000 is recommneded. 
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FIGURE A-IIIB 


SAMPLE OUTPUT 


[diagonal case] 


FIGURE ID¢o————398 .4000 Print out of Regis- 
2.0413 ters 0-9 and A-I 

O.2173 

-0.2590 

ball25 

-0.1971 

0.1904 

0.0027 

-0.0084 

0.2636 

-0.4521 

22 2 ea 

£ yf 13.5882 

10.0860 

1.0336 

[AW DN we 1.5332 
P ee 991.4751 


! 


——— ee 
RE, /RE 


HMoOnwruwuo@DysnuURWNFO 


0.8350 
3.1448 
0.0001 
0.0813 
-0.3264 
Seo 0 
Os Syeda: 
Onis 260 
-1.0000 
Vo290 


Print out of Regis- 
ters 10-19 


wo OsI A US WN FO 


Value of REw being 


evaluated Pe 10000. 0000 


3.743772534+10 oe / VOS 
PV og To 1. 370584009410 
320.2228 <———_—___—— On / Von 
ND A 

NTU /Vos 52.9601 
20000.0000 

2.502842147+11 

9.564330240+10 

50725 052 

80.4954 
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APPENDIX II 


The purpose of this appendix is to present: 


(1) A full description of the geometry of each sur- 


(2) 


a) 


face evaluated in this study. 


Friction factor (f) and Colburn Modulus (j) as 
a funtion of Reynolds number (RE) for each sur- 


face evaluated in this study. 


A complete compilation of the surface perfor- 
mance parameters required to compare each sur- 


face evaluated in this study. 
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APPENDIX III 


SIZING CROSSFLOW FINNED TUBULAR HEAT EXCHANGERS FOR A GIVEN JOB 


Crossflow finned tubular heat exchangers have found app- 
lication as gas turbine plant intercoolers, air conditioning 
unit heat exchangers, aircraft engine coolers, and numerous 
other uses. These applications most often involve gas-liquid 
service. 

The "job" of the heat exchanger will be defined as trans- 
fering a specified amount of heat between two fluids at given 
flow rates (W) and with specified amounts of pumping power on 


each side. The following values are specified: 


W dee ’ T 


i’ in, APy, ae : Y 


1 


Properties 


out, ’ ia 1 


1 


W ec oi ing a ’ AP 


2° in, out — : p 


5! in, oh! PYOperties 


2 


Subscript 1 will refer to the hotter fluid inside the 
smooth tubes and subscript 2 to the colder fluid flowing over 
finned tube banks. 

While core pressure drop accounts for by far the greatest 
percentage of the total pressure drop, entrance and exit loss- 
es would have to be included in the final design. [Ref 1 Chap- 
ter 4] 


Figure A-l shows the heat exchanger arrangement. 





SIDE 1 
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FIGURE A-1l 





~~~ Fin pees 
(1) In the air-conditioning industry, | 
fins which are common to all tubes in 
the bundle as distinct from having | 
each tube separately finned are Pree 
eee This is called plate finning. | 
| 
| 
| 


puter L-fins (3) Extruded (4) G-fins 
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The Basic Equations 


The first portion of the procedure involves determination 


of A. p A. » X, Y, and Z based on the specified values of AP, 
1 2 


and AP, and an assumed value of RE. - Then based on heat trans- 


fer considerations, it is determined if the proper heat balance 


exists. If not, another value of RE. is assumed and the process 


is repeated. 


4 Ae nt D 
= ie Be 
By = = [| D =D... (1) 
Pah AY Fob fo ound x a 
-T Sass 
4A 
. 
B = ———— = A specified geometric constant for (2) 
Doh cu Surface 2; See Figure A-II 
D5, = 4Fy! 
A 
c D 
mol 2 ih, (3) 
i 1 
Agr 4 
1 
A 
Cc D 
2.2, 2h. x (4) 
A 24 2 
fr. 
D 
aA = (6. —2B) [xz] = K, Xz (5) 
Cc) 1 4 1 





From 
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D 
A, = [By <7 [YZ] = K, Y 2 (6) 


| ee eee ee (8) 


equations (5) and (6): 


, (Wy IK, ¥2] 


G 
— Fe K, uf (9) 
G. [W.] ([K, XZ] Z 
WwW, K 
where K = 
ey 
2 
Din 7 Jo Om 


where L is the flow length: 


Side l: L= Y 


It 


Side 2: L X 





From 
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[AP, Dy, 0,3] [£5] [XI] 


- ee a (11) 
s 
Go «EAPS Dan Oma! ff) IYI 


[AP] Dan Omi! 


Define Ky = pe 
4 APS Don Oma! 


Equations (9) and (11): 








2 
G [£,) [xX] Z 
+ ax, ae KS ee (12) 
G5 [f,] [Y] [X] 
or 
2 3 
fo ; K, [Y] —_ oa (13) 
3 5 
f) K, [X] [X] 
K; 
where Ke = 
K 
4 
Eliminate [Y/X] with equations (9) and (13): 
3 3 
f. 4 [G,] [K,] — [G,] (1.4) 
3 3 616 13 
fy [Go] [K,] [ > 





3 3 
a , el a 
6 3 3 
1 [RED U9 
3 
> REQ 
RE. £. = 
m1 K 
6 
Solving for RE, : 
a RE? Te 
2 2" 2 
RE 
1 f. Xk Te 
6 i 


To begin the deSign procedure, a value of RE. 


Use the f. 


verses RE plot for side 2 to determine f.. 
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3 
D } 
= (15) 
D ] 
Lh 
ates 
u D 
a" ih mie 
ae 
Hy 2h 
43 0.333 
= (16a) 
D 
2h 


1s assumed. 


A sample 


data plot from reference [2] along with the associated geomet- 


ric parameters describing a typical side 2 flow arrangement is 


shown in Figure A-II. 





99 
ELGURE A=11 


eel ee ee nn 


sume 


ans7 | aoe 
Te Sees 





N, x 10° (47,G/p) 
2.0 30 


CF- 8.72 


FIG. 93 


VINWED CIRCULAR TULES 
SURYACE CP - 6.18 


Tube eataide diameter - 0.38 is. 

Fin piteh - 8,72 per inch 

Plow paesage bydraulic diameter — 47,00.01288 ft. 
Pin thickneee (averags)*® ~ 0.018 ia. 

Free-flow aren/frentel area ~.¢ of. 824 2,..3 

Heat tremefsr sree/tetal volume - H al43 £t,/ft. 

Fin area/tetal area - 6.010 


Motes Experinentel uncertaiaty for hsat transfer 
reeults poesibly somewhst greater than the 
soninal o5% quoted for the other surfsces 
because of the wmeceeaity ef sstinating o 
ceatact resietance is phe bl-aetal tabee. 


Bw Pine slightly tapersd. 
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It should be noted that the friction factor verses RE plot 
for flow over finned tube banks can be closely approximated by 
a straight line on a Log-Log plot. This implies that f. can be 


accurately represented by an equation of the form: 


B 


f., = A, [RE] 


ie 


Table I list the linearized equations of f., [and J5] for 


2 
Figures 92-101 of reference [2]. 
This simplification allows for ease in programming as dem- 


onstrated in the suggested design procedure flow diagram out- 


lined in Table II. 


For exact calculations, fy) should be read from the f, ver- 
ses RE plot for the specific tubular surface under investiga- 

tion, if such data is available. From this plot for surface l, 
it is possible to determine RE, and fy to satisfy equation (16) 


or (16a) once a value of RE. has been assumed. 

If, after assuming a value of RE,» the required value of 
fy and RE, lay outside the range of data for surface l, a new 
value of RE, must be assumed. 

If the required data for surface 1 is not available, fric- 


tion factor for turbulent flow inside smooth tubes is well cor- 


related and given by Nikuradse as: 


= = 4.0 Log [RE vf] - 0.40 (17) 
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The transition from laminar to turbulent flow inside smooth 
tubes Occurs between Reynolds number of 2300 and 4000. 

Excellent linear approximations of the Nikuradse relation 
are given by: 


= Oars > 
l 0.0791 RE 2300 < RE < 20,000 (18) 


rh 
2 


=U 7 
l 0.0460 RE 20,000 < RE < 10 O19") 


rh 
2 


In a crossflow finned tubular heat exchanger the dominant 
heat transfer resistance will probably be the gas outside the 
tubes flowing over the fins (side 2). It is for this reason 
that high accuracy for the £ (and j) prediction of the fluid on 
Side 1 is not critical. 

Note that equation (16a), uSing the linear approximation 


for fy given by equation (18) may be rewritten as: 


- 3636 
Sis ees 
EAYRE. uw. DBD 
RE, = 64 0 eee ee (20) 
3 .3 
0.0791 Ke Wy Don 
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Having satisfied either equation (16) using f, verses RE 


Ii 


1 


data for the specific tubular surface under investigation, OR 


equation (20) using the specified linear approximation for tur- 


bulent flow inside smooth tubes, length dimensions X and Y may 


be calculated from the relations given in equation (10): 


3 
ee ee: 
- 2 £. RES u? 
ae? ale 
AP. D> 0 
Mee 1 th -o "1 
— > 
2 fy RE} Wy 


Length Z may be calculated from equations (5) and (7): 


4 Wy 


RE Xx 


ooh can Beal 


The number of tubes can be determined from: 


W, D 





a oes le 
ie ha 
aT Wy A. 


1 


As was the case for friction factor, Jy should be read 


(22) 


(23) 


(24) 


(25) 


from the I verses RE plot for the specific tubular surface un- 


der investigation, if available. 
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In the event that the actual heat transfer data for side 
1 is not available, correlations can be used without introduc- 
ing significant deviations as explained earlier. 

The Colburn correlation for forced-convection, turbulent 
flow in tubes is widely accepted: 


=0\520 


j, = 0.023 RE (27) 


The heat transfer coefficient [{h] may be defined as: 


ee Cc ian RE 
h, = a (28) 
-667 D 


eer Lh 


Or uSing the Colburn correlation expressed by equation (27): 


0.80 
0.023 RE Guess 
h, = i... oe) (28a) 
pr 667 p 
1 Lh 


The value of J5 should be read from the J5 verses RE plot 
for surface 2. As was the case for f,, the Colburn modulus can 
be expressed in an equation of the form: 

= 
Jo = A; RE 


Table I presents the linearized equations of J5 for Fig- 


ures 92-101 of reference [2]. 
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5 u 
oo Se Ea: (29) 


Refering to Figure A-III, calculate the appropriate value of 


fin and surface efficiency for side 2: 


ee eee ee Bk; 
AU A U A U h A n h 
ie iy mie 1 eee, 22 
where: 
Ay, = a4 XYZ 
1 
An = a5 XYZ 
2 
If al Cor Cin me Cie Sonnac = C,/C 
"in, a Tout, 
(FS (34a) 
Bee ee 
1 2 
If Co Cys Cain - Cor Cain’ °mex a Co/Cy 
qT - T. 
out. in, 
PE ee eS (34b) 
T - 7T 
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FIGURE A-III 


Compute fin efficiency assuming a 


Longitudinal fin of Rectangular 


Profile: 
_ Tranh mé 
a me 
30 
ZOU ( ) 
where 2h 
2 
m 


What aeometric shape describes the 


fin being employed on the tube 


surface? 









Radial Fin of Rectangular 
Profile 


Tapered (Hyperbolic) Rad- 
ial fin 


Rectangular Profile fin 





sheet 





{see below for 2 defn] 


















An 


&n 0 ‘nN correction 


ea = 
correction correction 





Cy + Dy oO 


Rectangular Profile 
* 





whera: using r (staggered bank) 
.50 
ae [23 ead 
e T Z 


re) 
C,= -.00259 =- .0084 ln — 
4 ar 









r 
D.= ~.02426 - .4521 ln -2 
A ro 


[Ref 3 page 264] 











* 
L=r er 
e re) 


[Ref 4 pg.3-116] [Ref 4 pg.3-116] 


Vi-corrected — 


correction (31) 





Compute Surface n: 


Aq 


0 . "corrected! 
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also 





NTU = (35) 


If side 2 is in the form of rectangular fin sheets, the 
gas flowing over the tubes and the fluid in the tubes are con- 
Sidered UNMIXED. Using Figure A-IV with C_._/C and the mag- 

SE min’ “max 
nitude of € calculated from equation (34a or 34b) read the mag- 
nitude of NTU required to produce this effectiveness. Find 


[NTU Te and compare with the results of equation 


required! !"mi 
(33). If the equation is not balanced, a different value of 


RE. must be assumed and the procedure repeated. 


If side 2 is in the form of low fins, the gas flowing over 
the tubes is considered MIXED while the fluid in the tubes re- 
mains UNMIXED. The procedure outlined above is used with Figure 
A-V to find the NTU required to produce the calculated effect- 
iveness in this case. An alternative to using Figure A-V is to 
calculate the NTU requirement from equation (14) or (15) of ref- 


erence [2]. 


Table II presents the entire design procedure in a flow 
diagram. 
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FIGURE A-IV 


[BOTH FLUIDS UNMIXED] 








EFFECTIVENESS € 





NO. OF TRANSFER UNITS, NT Umox? AU/Gmin 
FIGURE A-V 


[ONE FLUID MIXED; ONE FLUID UNMIXED] 


EFFECTIVENESS 





NO. OF TRANSFER UNITS, NTUpmox= AU/Gmin 
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TABLE I 


B 


LINEARIZED EQUATIONS OF THE FORM A [RE] FOR FRICTION FACTOR & 


COLBURN MODULUS FOR FINNED TUBULAR SURFACES PRESENTED IN REFER- 


ENCE [2] 





COLBURN MODULUS [3] 





FRICTION FACTOR [f£] 









FIGURE A B A. B. 






£ f 





NOTE: These equations are written for the range: 
1000 < RE < upper limit of data IAW Ref [2] 








ant a 





Secepitar Pcegeg 


Tir JIL pLieiiaei 


a —- 
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CROSSFLOW FINNED TUBULAR HEAT EXCHANGER DESTGN PROCFRDURE 
ee en ns ern Rr SS eS i enero 


INPUT 
Linearized equation for f and j of 


the form A{RE}® for the enhanced 
surface: (subscript 2) 


Aecey? Bewey? Az(ey® 5 (e) 


Subscripts: 


Geometric Factors 1 Hotter fluid in 


Smooth Tubes 


2 Colder fluid 
flowing over fin 
ned tube banks 


) ae 
Dine Pane By Bo. Se ae 


Z, & 
Design Constraints: 
SP,. P 
¢ 





W T T 
ps in? out, in, 


W T 


2, T AP Pe 


in, out, eli in, 
Other: 
Properties of the two working 
fluids, fin material [{k] and 
shape, constants used to cor- 
rect fin efficiency 





bP. Po 8, We 


K.= po 
6 AP, O, Bo Wy 


Assume a Assumption: 
YD value for RE. 


0.25 


£,=0.0791 RE, 





34B 
Agia) RE. 
(.0791) Ke u 


f (e) 


RE,= 
i 
1 D 





TABLE II 
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Go To B to compute 


n 
6 


Calculate: Calculate: 


NO a Ce YES 
Go To A and Assume End of design Procedure 
another value for 
RE 


2 








lil 


Compute fin efficiency assuming a 
Longitudinal fin of Rectangular 
Profile: 


Tanh m2 
ur os me 


-50 


where: 2 
: se 


What geometric shape describes the 
fin being employed on the tube 
surface? 


Radial Fin of Rectangular Rectangular Profile fin 


Profile sheet 
(see below for 2 defn) 


An An An 


= = 0 Z = 
correction correction correction 


Rectangular Profile 
t 


whera: using r (staggered bank) 


50 
e T 2 


D,= --02426 - .4521 In : 
e ~ fo 


[Ref 3 page 264) {Ref 4 page 3-116) [Ref 4 page 3-116] 


"f-corrected e 


NeHL ‘Ncorrection 









Compute Surface n: 











] 


ry -corrected 


A,:; 
a = _ fin in - 
n 1 Re [ Lie 





iba ly 
TABLE III 


SUGGESTED PROGRAM FOR SIZING CROSSFLOW FINNED TUBULAR HEAT EX- 


CHANGERS FOR A GIVEN JOB IAW TABLE II* 


I. Design Constraints/Geometric Constants to be stored in 


Memory as follows: 


Primary Regis ces Secondary Register Other 
Quantity Units Quantity Units Quantity Units 


me /Bae BTU/LBM-F RE. (assumed) ND 


2 2 


ND Go?m2  LBM/HR'-FT 0.00 


ND Wy LBM/HR 0. 
LB/HR-FT FT? /Fr? 


Pa LB/FT? 


ND LBM/HR?-FT? 
ND 

LB/HR 

LB/HR-FT 


FT 





ioe. OUTPUT 


ee SIO Anko be applied’ "{-corrected’ Tene Cnin’cal'td 


*Por Hewlett-Packard HP-67/97 Programmable Calculator 





Uta) 


mol. PROGRAM FOR SIZING FINNED TUBULAR HEAT EXCHANG=RS FOR A GIVEN 





® 
JOB 


[Table III continued] 


857 
e658 
859 
668 
061 
G52 
063 
064 
@o5 
666 
067 
868 
869 
878 
67! 
872 
873 
074 
875 
076 
Q77 
#78 
879 
e890 
e861 
882 
683 
084 
885 
686 
087 
688 
e39 
898 
Q5! 
892 
893 
694 
895 
696 
897 
898 
899 
166 
161 
102 
163 
164 
165 
180 
187 
108 
169 
118 
fi 
bi2 


53 
-24 
oo fi 
vc 

3o bE 
-59 
3! 
24 
82 
~24 
36 a! 
~é4 
30 U2 
Os 

of 
#35 
16-51 
30 8! 


@S7 


2 On Te ms GN Gs 


GI 
q 
« Gy Oe Sa me CH MG 


LI 
nw 
4 


=e 


113 
1]}4 
Hee) 
116 
hiv 
113 
119 
128 
12! 
122 
123 
124 


147 
148 
149 
15¢ 
15! 
132 
153 
154 
155 
156 
137 
158 
153 
168 
Jol 
162 
1o3 
led 
165 
1eé 
lo? 


Lex 


Fes 
ReL4 


A 
RCLS 
x 


PRIA 


AcLe 

x 
$70C 
RCLA 
RULS 
! 

? 
y* 
Role 


*For Hewlett-Packard 67/97 Programmable Calculator 


IV. OPERATION 


First Run: 


Subsequent Runs: 


Select a Value of RE 


A. 


required value of NTU C_. 
min 


and 


1 
° 


3 
! 
3 

is 
5 


oO, 


3 
! 
3 
3 


oO. 





store in register 


To begin the program, press key "A" 


[As necessary until the calculated and 


value of RE 


begin the pfo 


agree] Select a 


and store in register A. To 
gram press key "B". 
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HEAT EXCHANGER DESIGN PROBLEM 


A. Description 


APPENDIX IV 


Appendix III procedure can be used for sizing crossflow 


finned tubular heat exchangers. 


The following data, taken from 


reference [2], is used to determine the required heat-exchanger 


Size for the given conditions. 


The design is that of an intercooler for a 5000 SHP gas 


turbine plant. 


DESIGN CONSTRAINT 


Surface 


SEDER SL 


Smooth Tube 


01224 FT 
N/A 

42.1 PT? /Ppr? 
N/A 
N/A 

400,000 LB/HR 
60 °F 
82 °F 

1.554 PSI 

2.36 LB/HR-FT 


1.0 BTU/LB-°F 


62.38 LB/FT” 


SIDE II 


Finned Tube 11.32- 
.737 SR [Fig 106} 


Otto c ate 
O0033 FE 
270 FT*/Fr° 
~845 
.225 IN 
200,667 LB/HR 

260 °F 

80 °F 

.306 PSI 
.0482 LB/HR-FT 


.2436 BTU/LB-°F 


-1755 LB/FT> 





i lis: 


SIDE I GIDE “Ik 
oT 6.8 .70 
K a N/A 100 BTU/HR-FT-°F 
F. .55 IN N/A 
b .79 IN N/A 
A -0560 IN? Calculate 


(From Fig 106) 


B. Calculations 


_ [.0560 IN‘] [4] 


A 
1 fro» [.55 IN X .79 IN] [.01224 FT] 


= 42.1 Pr? Pr? (1) 
B, = 270 Pr’ Fr? (2) 
K, = [42.1] [.01224/4] = .129 (3) 
K. = [270] [.01152/4] = .778 (4) 


qk. = £400,000] [.778] ~ 15.038 (9) 


(2007667) [.129] 


K. = [1.554] [.01224] [62.38] _ 1917.8 (11) 


PSvioele tt. Ollo2) 9 f.1755] 


=> 
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2 

Xe = 112.038)" _ 756 (13) 
[1917.8] 

K, = 1.0756] qo llaseey x 10"> (14) 
[12.038] 


CALCULATIONS FOR THE FINAL RE. SELECTED WILL BE SHOWN: 


Select RE. = 5760 


From the data of Surface 2 [Figure 106] read: 


Jo = .0054 

In order to be strictly consistant with the design presented in 
Ref [2], friction factor (f,] will be taken from Figure 29 of 
Ref {2], rather than employ any linear approximations suggested 


in the design procedure: 


From Figure 29 data for surface 1 and Equation (16): 


RE, = 5073 fy = .0072 
RE, wu 
G. = _ 2 2 = {5760} [.0482] = 24,100 LB/HR-FT“ 
Do, {.01152] 
2 
es {.306] [.01152} [144] [32.2] [3600] [.1755] (227 


[2] [0.021] [24,100] 7 





dys 


X = 18.29 IN (22) 
we 5.34 FT (23) 
a = [400,000] /978,127] _ ee (24) 


p29) 18.267 12) 


(5073] [2.36] 
[~01L224) 


Z 


where Gy = = 9, L2/ .B/HR-FT 


peertis.291 124.96) _ 55) (25) 


pO So 10. 79) 
Once again, in order to be strictly consistant with the design 
presented in Ref [2], Colburn Modulus [3,] will be taken from 
Figure 28 of Ref [2], rather than employ any linear approxima- 


tions suggested in the design procedure: 


From Figure 28 data for surface l @ RE, = 5073: 
Jy = 0.00502 (read Nee = .0014 and compute j] 
h, = [0.005] [1.0] [978,127] 16.8) 7°09! 
= 1368 BTU/HR rr? oF (28) 
-.667 
h, = [0.0054] [.2436] [24,100] [0.700] (29) 
= 40.21 BTU/HR PT? OF 
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m, = (2) (40.21) (.00033) (100) 
_ -] 
ee on Tr (30) 
Near mecanin | 295 se | 22572) 2 0.788 (30a) 
fare 3) (1.225712) 
mee eetion = 0 [For consistancy with Ref (2) (31) 
example] * 
ie = 1 - .845[1 - .788] = 0.82 (32) 
i ere 
er ealeuiated see )ii7 11.55] [40.21] [0.82] [4570.9] 
Z ee ee (33) 
130,506 BTU/HR-°F 
In order to complete the design procedure, compute AU equired 


and compare with equation (33) results: 


[AU] NTU C_, 


required in 


@: LJ = C 


min air (side II) [20076607 |) [12436 | 


il 


48,883 BTU/HR-°F 


E Z x: 
max = Cwater(side 1)~ 490,000 BTU/HR-'F 





Vay 


_ 260 - 82 


= 0.89 (34b) 
260-7 60 


From Figure 5 of Reference [2], for a cross-flow exchanger with 


both fluids UNMIXED, and C Vie = 0.122 with « = 0.89, for 


min’ “max 


from Table 4 of Ref [2] with proper interpolation] the eae 


is found to be 2.67. 


NTU = [2.67] [48,883] = 130,517 BTU/HR-°F (35) 


req' amin 
Comparing Equation (33) with equation (35): 


5 
NTU = NT 


req'd°min Vealculated Cin 


130,517 BTU/HR-°F 


130,506 BTU/HR-°F 


CLOSE ENOUGH! 


Surface 2 1s a rectangular fin sheet and, to be strictly 
correct, a correction should be applied to permit the use 
of equation (30a) as outlined in the design procedure. [See 
Appendix III, Table I] In order that the results of this 
example will be consistant with those of reference [2], a 
correction of 0 will be assumed pursuant to the simplifying 


assumptions as stated in the example of reference [2]. 
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